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We demonstrate a significant resolution enhancement beyond the conventional limit in multipho¬ 
ton microscopy (MPM) using saturated excitation of fluorescence. Our technique achieves super- 
resolved imaging by temporally modulating the excitation laser-intensity and demodulating the 
higher harmonics from the saturated fluorescence signal. The improvement of the lateral and axial 
resolutions is measured on a sample of fluorescent microspheres. While the third harmonic already 
provides an enhanced resolution, we show that a further improvement can be obtained with an 
appropriate linear combination of the demodulated harmonics. Finally, we present in vitro imaging 
of fluorescent microspheres incorporated in HeLa cells to show that this technique performs well in 
biological samples. 


INTRODUCTION 

Since its introduction in 1990 [l[, multiphoton 

microscopy (MPM) rapidly evolved to become the 
workhorse for 3D fluorescence imaging in many labora¬ 
tories. MPM is arguably one of the most influential re¬ 
cent inventions in optical microscopy because it offers a 
combination of minimally invasive imaging with high res¬ 
olution, thin optical sectioning and deep penetration in 
highly scattering samples such as biological tissues (H . 
These properties all stem from the process of multipho¬ 
ton fluorescence excitation. In order to achieve efficient 
multiphoton excitation, photons need to be concentrated 
both in space by focussing the excitation beam, and in 
time by using pulsed laser sources. Fluorescence is there¬ 
fore only excited in a narrow volume. Moreover, because 
the excitation spectra of common fluorescent molecules 
are in the visible spectrum, the wavelengths used for 
MPM are in the near-infrared such that optical sample 
damage is much reduced [1, . 

A second important development in optical microscopy 
has been the development of super-resolution microscopy 
which achieves resolutions beyond the diffraction limit. 
In essence, super-resolved imaging requires the sample 
to respond non-linearly to the illumination light. Such 
a non-linearity can be achieved by, among other tech¬ 
niques, saturation of excitation Q and stimulation emis¬ 
sion depletion (STED) Bringing super-resolution to 
MPM has the potential to enable studying biological 
samples at the depths offered by MPM with a resolution 
beyond the diffraction limit. 

STED has been applied successfully for two-photon mi- 
crosc^y (2PM), but only to improve the lateral resolu¬ 
tion [Tj. Moreover, STED requires the use of a high- 
intensity visible wavelength depletion beam, increasing 
the invasiveness of the method. Although STED mi¬ 


croscopy can also improve the axial resolution beyond 
the diffraction limit [8|, it requires a carefully aligned 
and relatively complex optical system, often including 
an additional depletion laser source. 

In this work, we present a simple approach to achieve 
3D super-resolution in MPM, based on saturation of flu¬ 
orescence excitation Q. This method can easily be im¬ 
plemented on an existing multiphoton microscope, as it 
only requires (a) temporal modulation of the excitation 
beam using an off-the-shelf intensity modulator, and (b) 
demodulation of the detected fluorescence signal with an 
electronic lock-in amplifier. Hence, no additional laser 
sources or beam shaping optics are required, and only 
near-infrared light is incident on the sample. 

Images of fluorescent microspheres on a coverslip prove 
that this technique can be implemented successfully in 
MPM and improves not only the lateral but also the ax¬ 
ial resolution. Eurther experiments with microspheres 
integrated in HeLa cells produced encouraging results 
demonstrating the feasibility of this method for in vitro 
studies. 


PRINCIPLE OF THE METHOD 

The technique proposed in this work is based on the 
saturation of the probability of excitation of the fluo- 
rophore in MPM. Until recently, saturation has been 
considered harmful for imaging and has been avoided 
It is possible, however, to show that in a satura¬ 
tion regime, the relation between the intensity of excita¬ 
tion and the emitted fluorescence intensity is no longer 
linear and contains nonlinearities that can be used to im¬ 
prove the lateral as well as the axial resolution iinilii- 
In a two-photon excitation process, assuming a delay 
between pulses much longer than the fluorescence life- 
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time, the probability of exciting a fluorophore at a given 
position with one excitation laser pulse is given by 

= 1 — ^-(^cr2PAP^PSF^{x,y,z) 


where a is a coefficient which depends on the pulse’s 
shape and duration, and on the repetition rate of the 
laser, (J 2 pa is the two-photon absorption cross-section of 
the fluorophore, P the average power of the excitation 
beam and PSF{x,y, z) is the illumination point spread 
function. To extract the nonlinearities that appear when 
saturation occurs, a temporal modulation of frequency 
is applied to the excitation beam @ . The probability 
of excitation can then be written 

p _ 2 _ ^-Oicr2PAP‘^[l-\-cos{ujt)]‘^PSF‘^{x,y,z) / 2 ^ 


where uj = 

At high excitation power, higher order harmonics in 
the fluorescence signal will become non-negligible, as il¬ 
lustrated in Fig. [IJa). 

By developing the 1 — e~^ term up to the 4^^ order in 
Eq. (|2]), and rearranging the expression, the power in the 
harmonics can be calculated 

( p^^anai . 2PSF‘^ - iPSF^ + ^PSF^, 

< ■ IPSP"^ -iPSF^ + ^PSF^, (3) 

. PtT' ■ \PSF^ - §PSF^. 

These equations show that the harmonics in the de¬ 
tected fluorescence signal contain terms with the fourth 
and sixth powers of the illumination PSF. In stan¬ 
dard two-photon microscopy, the detected signal contains 
PSF^. Hence, the resolution is improved by at least 
a factor of \/2 when using the third harmonic. Figure 
[TJb) shows the dependence of the harmonics’ strength 
Semission OH the illumination intensity. 




Frequency (kHz) log(a2PAP^) 


FIG. 1. Apparition of higher harmonics in the saturation 
regime, (a) Spectrum of the signal of the detected fluorescence 
at a modulation frequency of I'm = 5 kHz. (b) Power in the 
n-th harmonic of the detected fluorescence signal (Semission) 
for increasing power on the sample. 

By making the appropriate linear combination of the 
harmonics in Eq. it is possible to obtain a signal 
which is dependent only on the sixth and higher powers 


of the PSF: 


Linear combination signal = p^^Qnai_^ps^gnaips^gnai^ 

^4) 

The linear combination signal (oc PSF^) therefore 
achieves an additional resolution improvement by a factor 
of y/3j2 over the third harmonic signal assuming 

a gaussian illumination profile. 


MATERIAL AND METHODS 


Experimental setup 



Processing 


FIG. 2. The experimental setup is based on a conventional 
multi-photon microscope with an acousto-optic modulator to 
modulate the fluorescence excitation intensity, and a lock-in 
amplifier to extract the harmonics of the fluorescence signal 
detected by the PMT. 


The experimental setup, shown in Fig. [2l is built 
around a commercial two-photon microscope (Model 
FVIOOO-MPE, based on a BX61WI upright microscope 
body, Olympus Corp.). The configuration of the micro¬ 
scope is the U-MWG2 which has two detection chan¬ 
nels with filters allowing the detection of fluorescence 
from 510 to 550 nm and from 600 to 690 nm. The mi¬ 
croscope objective is the Olympus XLPN25XWMP (nu¬ 
merical aperture 1.05, magnification 25x). The multi¬ 
photon excitation laser system used in our experiments 
is a Spectra-Physics Ti:Sa Tsunami pumped by a Mil- 
lenia ProS-Series. This laser is tunable over a wave¬ 
length range of 840 nm to 980 nm. It has a repetition 
rate of 82 MHz and is capable of producing 150 fs pulses 
at a wavelength of 880 nm. The excitation light beam 
is temporally modulated by an acousto-optic modula¬ 
tor (MT80-A1-IR, A A Opto Electronic Sa.). A lock-in 
amplifier (HF2LI, Zurich Instruments AG) demodulates 
the signal from the photomultiplier tubes (R3896, Hama¬ 
matsu Photonics K.K.) of the microscope at the first, 
second and third harmonics. The linear combination of 
harmonic signals is generated by an analog summation 
circuit such that no offline processing is required to gen¬ 
erate the final image. The microspheres used for the 
different experiments are Dragon Green fluorescent mi¬ 
crospheres (Bangs Laboratories, (480, 520)). 
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Sample preparation 

Stock cultures of Hela cells were maintained in RPMI 
(Roswell Park Memorial Institute medium) supple¬ 
mented with 10% (v/v) deactivated fetal bovine serum, 
2mM L-glutamine and 50/ig/ml gentamycin (complete 
medium, Lonza) at 37°C in a humidified atmosphere of 
95% air/5% CO 2 . Cells were recovered by trypsinisation 
(trypsin-EDTA, Lonza). 

Fluorescent microspheres (Bangs Laboratories) were 
used, with a mean diameter of 0.51 /im (FC03F). The sus¬ 
pensions of beads were homogenized (rotation) overnight 
at 4°C before use. 5 x 10^ cells were mixed with 10 jjI 
of bead suspension and 4 ml complete culture medium, 
plated in 10 cm^ culture wells and incubated for 24 hours 
at 37°C in a humidified atmosphere of 95% air/5% CO 2 . 
At the end of the incubation the cells were recovered by 
trypsinisation. To eliminate the beads not adsorbed by 
the cells, cells were diluted in 2 ml medium, put in a 3 /im 
pore size cell culture insert (Greiner Bio-one) and briefly 
centrifuged at 200 g. Cells were recovered by washing the 
insert with culture medium, plated on glass coverslips at 
a density of 5 x 10^ cells/cm^, and incubated overnight 
at 37°C under an atmosphere of 5% CO 2 . The coverslips 
were then rinsed 3 times with phosphate buffered saline 
(PBS), fixed with 4% Paraformaldehyde (PFA) during 
10 minutes and rinsed again with PBS. 


Acquisition parameters 

All imaging experiments were performed at a mod¬ 
ulation frequency Um of 17 kHz, and the first, second 
(34 kHz) and third (51kHz) harmonics were demodu¬ 
lated. The images containing only the first and second 
harmonics were obtained with 10 mW excitation power 
on the sample. To obtain the images of the third har¬ 
monic as well as the linear combination, the excitation 
power was increased to 25 mW to produce a sufficiently 
strong signal. The complete 2D held of view (x, y) of 
508.68 X 508.68 yw? is sampled at 4096 x 4096 pixels with 
a pixel dwell time of 200 /is, which corresponds to 3.4 
modulation periods at 17 kHz, 6.8 at 34 kHz and 10.2 at 
51kHz. {x^y^z) scans are stacks where each (x,^)-slice 
(same held of view as for 2D images) is sampled at 
2048 X 2048 pixels. 


RESULTS 

Microspheres on a coverslip 

To demonstrate the resolution improvement of our 
method, 0.19/im microspheres were imaged. First we 
determined the power needed to obtain saturation by 
measuring the harmonics in the fluorescence signal. To 


this end, we analyzed the signal from the PMT with 
an electronic spectrum analyzer (ESA) (8564E, Hewlett 
Packard), and increased the illumination power until the 
harmonics rose above the noise floor of the ESA as shown 
in Fig.[3fa). Note that, due to the narrowband detection 
of the lock-in amplifier, the actual noise floor of the imag¬ 
ing system is expected to be lower than the noise floor 
measured by the ESA (—70dBm). 



FIG. 3. Harmonics in the detected fluorescence signal from 
0.19/im microspheres on a coverslip. (a) Frequency spectra 
of the fluorescence signal emitted by the microspheres for in¬ 
creasing average power, (b) Spectral intensities of the three 
first harmonics in the fluorescence signal. Linear fits of the 
data are shown as black lines. The particle concentration was 
lOmg/mL and the frequency of the excitation modulation Vm 
was 17 kHz. 

The third harmonic clearly emerged from the noise for 
an illumination power of approximately 30 mW on the 
sample. At this illumination power, the signal-to-noise 
ratio was measured to be 15 dB and was sufficient to ac¬ 
quire an image. The demodulated fluorescence intensity 
at a frequency of Vrn (17kHz), 2z/^ (34kHz) and 3z/^ 
(51 kHz) is shown in Fig. ^h). As expected from Eq. (j3j) 
and Fig. [TJb), the curves of the fundamental and second 
harmonics are linear and have the same slope, while the 
slope of the third harmonic curve is twice as steep. 

For the fundamental (17 kHz) and the second harmonic 
(34 kHz) [Fig. H a) and (b)], the spatial resolution in x 
and y was measured to be the same as in conventional 
MPM. The images obtained using the third harmonic 
[Fig. |4l[c)] and the linear combination of the three first 
harmonics given by Eq. dH [Fig. md)] clearly show a 
resolution improvement. The intensity profiles shown in 
Fig.|4l[e) confirm this improvement: for the first and the 
second harmonics, the two microspheres are not resolved, 
they can just be resolved by the third harmonic. The 
linear combination resolves the two microspheres even 
better. 

The full-width at half-maximum (FWHM) was mea¬ 
sured on the leftmost microsphere to be approximately 
470 nm for the two first harmonics, 340 nm for the third 
and 285 nm for the linear combination signal. These mea¬ 
surements fit with the estimated theoretical improvement 
of the resolution. 

The same analysis was performed in the axial direc¬ 
tion to demonstrate the improvement of the axial res- 
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FIG. 4. Characterization of the lateral resolution improve¬ 
ment using 0.19 /im microspheres on a coverslip. (a-d) Images 
in the {x^y) plane reconstructed with the first, second, and 
third harmonics and the linear combination of these three har¬ 
monics. The (a:,y)-scan size is 48.36 x 56.54/xm^ (390 x 456 
pixels) resulting in a total acquisition time of ~ 35 s. (e) Sig¬ 
nal profiles along the orange line in (a-d). The background 
intensity was measured in an empty area and subtracted ev¬ 
erywhere. 


olution. Figure [5] shows the four different (x, z) planes 
obtained with the first, second, and third harmonics and 
the linear combination. The results confirm that the ax¬ 
ial resolution is also improved using saturation of the 
excitation. The axial FWHM was measured to be ap¬ 
proximately 3.4/xm for the two first harmonics, 2.1/xm 
for the third and 1.8/xm for the linear combination sig¬ 
nal. Again, these values match the estimated theoretical 
resolution improvement. 
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FIG. 5. Characterization of the axial resolution improvement 
using 0.19/xm microspheres on a coverslip. (a-d) Images in 
the (//, z) plane reconstructed with the first, second, and third 
harmonics and the linear combination of these three harmon¬ 
ics. The (//, z)-sections are taken from (x,//, z)-stacks consist¬ 
ing of 19 (a:,//)-slices of 49.10 x 58.53 /xm^ (192 x 236 pixels), 
resulting in a total acquisition time of 3 minutes per 3D stack, 
(e) Signal profiles along the yellow line in (a-d). The back¬ 
ground intensity was measured in an empty area and sub¬ 
tracted everywhere. 



FIG. 6. Microspheres of 0.51 /xm in HeLa cells, (a-d) Images 
in the (//, z) plane show a clear improvement of the resolution 
in the lateral and axial dimensions, (e-h) 3D renderings for 
the first, second and third harmonics and the linear combina¬ 
tion of the three first harmonics. The (x, //, z)-stacks consisted 
of 17 (a:,//)-slices of 169.6 x 169.6/xm^ (800 x 800 pixels), re¬ 
sulting in a total acquisition time of 36 minutes per 3D stack. 


Microspheres in HeLa cells 

To demonstrate the applicability of this multiphoton 
super-resolution technique to biological samples, we went 
on to image 0.51 /xm microspheres incorporated in HeLa 
cells. Figure [6] shows {y^z) sections of 3D images con¬ 
structed from the first, second, and third harmonics and 
the linear combination of those three harmonics. 3D ren¬ 
derings were generated with the Imaris software (Bit- 
plane AG). 

As for the microspheres on a coverslip, an improvement 
of the resolution in all three dimensions is observed when 
imaging at the third harmonic and a further improvement 
is obtained by the linear combination of the three first 
harmonics. 


DISCUSSION AND CONCLUSIONS 

In this work, we demonstrated experimentally that sat¬ 
uration of fluorescence excitation can be exploited to ob¬ 
tain a 3D resolution improvement in MPM. Implemen¬ 
tation of this technique on an existing MPM system is 
straightforward since it only requires modulating the ex¬ 
citation beam and demodulating the detected fluores¬ 
cence signal. The method is also cost-effective because 
it does not require additional laser sources or microscope 
objectives to improve the 3D resolution. Moreover, the 
technique does not employ delicate 3D beam-shaping op¬ 
tics, which are sensitive to misalignment over time. 

Images of fluorescent microspheres on a coverslip and 
integrated in HeLa cells show that resolution is indeed 
improved in all three spatial dimensions by a factor of 
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\/2 when imaging at the demodulated third harmonic of 
the detected fluorescence signal. Furthermore, we showed 
that the 3D resolution can be improved by an additional 
factor of y^ 3/2 when the appropriate linear combination 
is made of the first three harmonics. This linear com¬ 
bination of harmonics is made by an analog circuit such 
that no subsequent signal processing is required. For 
comparison, our super-resolution MPM obtains the same 
lateral resolution as a 400 nm illuminated linear confo- 
cal microscope (assuming the same beam characteristics 
and objective), while conserving the advantages of clas¬ 
sical MPM (i.e. near-infrared illumination and increased 
depth of penetration). 

To obtain the required saturation regime and get a 
sufficiently high signal-to-noise ratio for the harmonic 
signals, an average illumination power of 25 mW at 
880 nm was delivered on the sample. Interestingly, this 
rather modest illumination power suffices to obtain a 
near twofold 3D resolution enhancement in MPM. More¬ 
over, since the sample is illuminated only at near-infrared 
wavelengths, our approach is likely to be less invasive, in 
terms of photo-toxicity, than the super-resolution tech¬ 
niques based on STED. 

With the present acquisition parameters, relatively few 
modulation periods are measured at every image voxel. 
Even though the resulting signal-to-noise ratio was suf¬ 
ficient for imaging, there is room for improvement by 
increasing the number of periods of modulation of the 
excitation signal per pixel. This can be done by in¬ 
creasing the pixel dwell time (if the acquisition time is 
not critical) or by using a higher modulation frequency 
(but keeping a sufficient number of pulses per modulation 
period). An improved signal-to-noise ratio could allow 
using a lower illumination power while still obtaining a 
sufficiently strong signal to reconstruct the image. This 
could be particularly interesting to investigate for future 
work, especially for in vivo applications where the power 
damage threshold is critical. 

Theoretically, the resolution improvement is unlimited 
since higher harmonics than the third one can be used. 
However, as for all super-resolution methods based on 
saturation of the excitation or STED, the attainable res¬ 
olution is ultimately limited by the damage threshold of 
the sample. At high illumination powers, photobleach- 
ing and phototoxicity become critical in biological sam¬ 
ples. This problem can be reduced using fluorophores 
with high two-photon absorption cross section or quan¬ 
tum dots, provided that their fluorescence lifetime is suf¬ 
ficiently short compared to the modulation period. 
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